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Abstract 

A general formalism is used to express the long-range potential energies in inverse powers of the 
separation distance between two like atomic or molecular systems with P symmetries. The long- 
range molecular interaction coefficients are calculated for the molecular symmetries A, 11, and 
S, arising from the following interactions: He(2 ip)-He(2 ^P), He(2 ip)-He(2 ^P), and He(2 3p)- 
He(2^P). The electric quadrupole-quadrupole term, C5, the van der Waals (dispersion) term Cg, 
and higher-order terms, Cg, and Cio, are calculated ah initio using accurate variational wave 
functions in Hylleraas coordinates with finite nuclear mass effects. A comparison is made with 
previously published results where available. 
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I. INTRODUCTION 



Accurate description of the interactions between two excited atoms (or molecules) at long- 



range is fundamentally irnportant for studies of molecular excited state spectroscopy 



associative ionization 



iiport 

m 



and other coUisional processes 



5|, and is at the heart of 



several schemes for quantum computation [6|, 3]- At sufficiently large separations, the mutual 
electrostatic interaction energy between the two excited atoms can be accurately described 
using an expansion of the potential energy in inverse powers of the separation distance R. 
The terms describe the electric quadrupole-quadrupole interaction at order and the 
instantaneous dipole-dipole (e.g. dispersion) interaction at order 8| and higher order 
instantaneous multipole-multipole interactions at orders R~^ and 

Long-range interactions involving few-electron atoms are the only interactions that 
presently can be rigorously calculated with high accuracy. Different levels of approxima- 
tion are needed for the calculations of long-range forces for alkali-metal and alkaline-earth 
J. fl y . Sizeable discrepances be.wee. vanous ca.cu.a.ions i. *e l.e.a.u.e ean 



atoms. 



occur, as illustrated in the comparisons of Cq coefficients, for example, given by Zhang et 
al. 12|] for Li (2p)-Li (2p) and given by Yurova [13] for Na(3p)-Na(3p). For helium, it is 



possible to perform a highly-accurate ah initio calculation of atomic properties and long- 
range interaction coefficients. Such results could become benchmarks for eventual ah initio 
calculations of alkaline-earth atomic interactions. Alkaline-earth and other two-electron ex- 
cited P atoms are currently being studied as the optimal candidates for frequency-based 
standards and optical clock experiments jl^ . 

We had previously studied the long-range interaction coefficients C„ (with n < 10) for all 
He(n^5)-He(n'^'5) and He(n ^5)-He(n' ^'P) systems of the energetically lowest five states: 
He(li5), YLe(/2lS\Y{e{2^S), He(23p) and YLe{2^P) and the finite nuclear mass effects for like 



isotopes 



15 



16 



171] . In this work, we present results for more complicated set of interactions 



between two like isotope helium atoms with P symmetries. Degenerate perturbation theory 
is needed to derive the interaction terms for some of the terms. Section II introduces a 
general formalism for calculating dispersion coefficients between two like atomic or molecular 
systems of P symmetry. Section III presents numerical results of dispersion coefficients C5, 
Ce, C's, and Cio for the following three systems He(2 ip)-He(2 ^P), He(2 ip)-He(2 ^P), and 
He(2 3p)-He(2 3p). 
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II. FORMULATION 



In this work, atomic units are used throughout. At large distances R between two atoms 
a and 6, the Coulomb interaction treated as a perturbation to the two isolated atoms, 
is 



=0 L=0 



where 

VeL = M-in^,Lr'^'j:KiL TjfKcT) T^'Jip) . (2) 
In the above, Tjp{(T) and T^^J{p) are the atomic multipole tensor operators defined by 

T^'\a) = j:Q^^lY,,{a,), (3) 

i 

and 

^?(p) = E^.pt^i-.(p.)> (4) 



■th 



and Qi and <Tj are the charge and the position vector of the i particle in atom a, respectively. 
Similarly, qj and p are for the j*'* particle in atom b. The coefficient K^j^ in Eq. ([2]) is 



' + fiJ\L + fi. 



(5) 



and {i,L,---) = (2£+ 1)(2L+ I)---. 

Since the Coulomb interaction V is cylindrically symmetric about the molecular axis R 
or z axis 0, [l^, the projection of the total angular momentum of the combined system 
a-b along the z axis (with magnetic quantum number M), is conserved. Therefore, states 
with M = ±2, ±1, and are not mixed with each other, corresponding to the A, 11, and E 
molecular states, respectively. The A and 11 states are degenerate with respect to the sign 
of M and the degeneracy can not be removed physically in the free combined system a-b. 
Therefore, we only study the states with positive M in this work. 



A state 



For two like isotope atoms a and 6 in P symmetry, the zeroth-order wave function for 
the A state of the combined system a-b can be written in the form: 

vl/W(A,/?) = ^[vl>„jM„;o-)vl>„,(M,;p)+/3v^„jM„;p)vl/„^(M,;o-)], (6) 



where Ma = Mb = 1 are the magnetic quantum numbers, a is the normahzation factor, and 
P describes the symmetry due to the exchange of two initial states \l/„^ and \l/nj^. If two 
atoms are both in the same P state, then a = \/2 and /? = 0: if they are in different P 

n 

states, then a = 1 and (3 = ±1 |10| . 

According to the perturbation theory, the first-order energy is 



yW = (^(o)(A,/5)|y|^(°)(A,/3)) 



(7) 



where, after some angular momentum algebra, one gets (see also Ref. 



C,{A,P) = -A, + PA2, (8) 

Air 

x{^M\\j:<i^p'Mp,mn,{p)) , (9) 

A'jr ^ . 

A, = -|(vi/„„(^)||EQ«^'^2(^.)ll^n,(^))r- (10) 

The second-order energy is 

= ^ y: \('^^'^(^)\^\x-siLsM,;a)ur.,iL,Mup))\' ^ ^^^^ 

nsUt LsMs LtMt ^n^nt ~ Enjnt 

where XusiLgMs; (T)unt{LtMt] p) is an allowed intermediate state with the energy eigenvalue 
En^nt = + En^, and the prime in the summation indicates that the terms with En^m = 
En]ni, should be excluded. Substituting Eqs. ([1]) and ([6]) into Eq. ffTTj) . we obtain 

v'-'--i:'E i: /'^''fo, ■ (12) 

nsm LsMs LtMt ^ngnt ^riaUf, 

with 

Si = \{^>^SMa]<T)^>^SMb-,p)\V\xnXLsMs-,(T)unt{LtMt-p))\\ (13) 

B2 = {<^nSMa\(T)mnSMb]p)\V\XnXLsMs](T)uJnXLtMup)) 

X {^>nSMa\p)'^nAMb]<T)\V\XnXLsMs](T)Unt{LtMt,p)) . (14) 

After applying the Wigner-Eckart theorem, we have 

v^' Bi ^ ^ Di(£, L, f , L') . 

2^ 2^ 2^ -r, 77(0) ~ D£+L+^/+L'+2 ' \^^) 

nsnt LsM^ LtMt -t^nsnt ~ -t^nani, U'LL' 



with 

Di(£,L,f,L') = ^ E Gi(£,L,f,L',L„L,)Fi(£,L,r,L',L„Li). (16) 

LsLt 

In Eq. f fT6|) . Gi is the angular-momentum part and Fi is the oscillator strength part. Their 
expressions are 

La £ Lg \ 



-Ma Ma-fl, 



and 



/I L Lt \/ 1 e W 1 L' \ 



F,{e,L,e',L',Ls,Lt) = — 

^ 9ns;nana{'^si ^1 £ )9nt\ni,ni,{Lf, 1, 1, L, L ) (18) 



with 



(^,^0 ^(2Li + l)(2L2 + l) 
X (Li; a) 1 1 E g.a,^r,(^.) I IXn. (1^^; ^)) 

i 

(vl>„^ (L2; ^) 1 1 E ^?^^f (^0 1 IXn. {Ls-, <t)) , (19) 



X 



and ^En^na = -^n^ ~ -E'na? ^tc. For the special case where the two initial states \E'„^ and 
"^nt are the same and i = i', gns\nana reduces to the absolute value of the 2^-pole oscillator 
strength 



Similarly, we have 



y^' B2 _ -D2(^, L,i', L') 

nsut LsMsLtMt ^risUt ~ -^naUb U'LL' ^ 



D2{i, L, i', L') = ^ E G2{i, L, L', L„ Lt)F,{i, L, L', L„ L^) , (22) 

LsLt 

with 

G2(£,L,f,L',L„L,) = (-l)^+^'(£,L,f,LO^/^E^^L^£'L''*'""'^ 
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4' ' ](' ' ]. (23) 

\-Mb -II Mb + iiJ\-Mb Mb -Ma + 11 Ma -II J 

F2{i, L, L', La, Lb, Lg, Lt) = 

2 n,nt (A£'„^^„^ + AEnt,nb)\/\^En^,na^En„nb^En^^na^En 



^nt,nb I 



Finally, the second-order energy V^"^^ is 



^-^--E^^^. (28) 



n>3 

where the dispersion coefficients C2n(A,/3) are defined by 

C2n{A,P) = Yl [D,{e,L,e',L')+pD2{e,L,i',L')]. (26) 

e,L,e',L'>i 

e+L+e'+L'+2=2n 

B. n state 

For the H state, the zeroth-order wave function is in the following form 

*W(n,/3,7) = -[^nSMa,CT)^nbiMb;p)+^^nMb;Cr)^nb{Ma;p)] 

a3 

+ ^[*„^(M„; p)^nb{Mb; (t) + j^nMb', p)*n,(M„; o")] , (27) 

where /? = ±1, Ma=Q, and M^ = 1. If Ua — rib, a. — a/S and 7 = 0. If na 7^ rib, then o; = 1 
and 7 = ±1. The first-order energy correction yields 

C5(n,/3,7) = -a2(l + 27 + 72)C5(A,/?). (28) 

The second-order energy correction is 

y(^)(n,/3,7) = -E' E E (29) 

Introducing a function Qo defined by 

g,{Ma,Mb,M„Ma,i,L,L„Lt,e,L') = {-lf+'^\l,L,e,L'Y'^ E ^^lK^v 

/I £ W 1 L LA/ 1 W 1 L' LA 

^\-Ma II Ms)\-Mb -II Mt)\-Mc fi' Ms)[-Md -fx' Mj' 
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one can write y(^)(n, /?, 7) as 



n>3 



(31) 



where 



C2„(n,/3,7)= Yl Dp,{e,L,i',L'), 

e,L,e',L'>i 

e+L+e'+L'+2=2n 



(32) 



D^^(^, L, L') = — E /i7(^' ^> Ci + /27(^, L, L„ Li, f , L') G2 

LsLt 



,(33) 



h^{i,L,L,,Lt,i',L')^ 



1 



(1 + 72)[^„(0,1,0,1,£,L, f,L') + (?o(l,0,l,0,£, L, Ls,Lt,i',L')] 



+jgo{0,l,l,0,e,L,L,,Lt,i',L')+jgo{l,0,0,l,e,L,L,,Lt,e',L' 



(34) 



/2^(£,L,L„Lt,£',L') = 

^(1 + 7^) [(7,(0, 1,1,0, L, Ls, Lt, L') + ^«(1, 0, 0, 1, L, L„ L*, f , L')] 
+7 ^„(0, 1, 0, 1, ^, L, L„ Li, L') + 7 ^„(1, 0,1,0, £, L, L„ L*, L') , 



(35) 



and 



UsTlt 



gns;nana{Ls, 1,1, i, i') gnt-mnjiLt. 1, 1, L, L') 

(AL;„,,„„+AK„nJ|A^„,,„„AL;, 



rit,n.6 I 



9ns\ni,ni,{Ls, 1, 1, i, ^')gnt\nana{Lt, 1, 1, L, L') 



(al;„,,„„ + A£„„„j|a^„, 



"6 I 



G2 = /5E' 



nsUt 



gns^riantiLs, 1, 1, I, f) gnt]ntna{Lt, 1, 1, L, L') 



gns;ntna{Ls, 1, l,i, i')gnt;nant{Lt, 1, 1, L, L') 



(36) 



(37) 

Prom the above expressions, it is clear that /i-y = 7/27 when 7 = ±1 and G2 — PGi when 
7 = 0. 
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C. S state 



For the S state, the possible zeroth-order wave functions for the combined system a-b are 
¥'\j:,/3) = ^[vI/„^(M,;o-)vI/„,(M,;p) + /?vl/„^(M,;p)vI/„^(M,;o-)], (38) 

and 

P, 7) = ^[^nSMa, o-)vl/„jM,; p) + 7*n.(M,; o-)vI/„jM,; p)] 

ad 

+ ^[vl/„jM,; p)^n,{M,; a) + 7^n.(M,; p)vl>„jM„; o")] . (39) 

In the ^(°)(S, /5) state we have Ma = Mh = 0, while in the ^(°)(S, 7) state, = -M;, = 
1, and 7 = ±1. When Ua = Ub, P = and a = -\/2, while otherwise P = ±1 and a = 1. For 
the S state, there is an additional symmetry with respect to the reflection of wave function 
on a plane containing the molecular axis. If Q is the reflection operator |lO|. 

g^,,(M„; r) = (-l)^^»vl>„^(-M,; r) , (40) 

then 

Q^(°)(E,/9) = ^(°)(E,/3), (41) 

and 

gxl>(°)(S,/?,7) = 7Vl/(o)(s,/3,7). (42) 

Due to this symmetry, the state \E'^°''(S, /3, —1) can not be mixed with other two states by 
the Coulomb interaction V. However, the degenerate perturbation theory is required to 
remove the degeneracy between ^W(S,/3) and ^W(E,/3,1). 

1. The first- order energy 

For the \E''^''^(S, —1) state, the first-order energy is 

(vI/W(S,/9,-l)|\/|M/W(S,/3,-l)) = 0. 

(43) 
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For ^'('^)(E,/3) and ^'(°^(E,/3, 1), the matrix elements of the Coulomb interaction V are 

{¥'\E,P)\V\^('\E,P)) = (44) 



and 



(vI/(o)(E,/3, l)|y|^(°)(E,/3, 1)) = -^^^^ . 



(45) 



i?5 

The first-order energies are obtained by the diagonalization of the following matrix: 

4C5(A,/?) 



(46) 



They are 



i?5 



Ai = 0, 

A2 = 



1 

1 1 

V2 2 



(47) 



6C5(A,/3) 



(48) 
(49) 



Their corresponding normalized eigenvectors are 



*(°)(E,/3; AO = -i=M/W(E,/3) + y|vl;(o)(E, /?, 1) , 
*W(E,/3; A2) = y|*^°)(S,/3) + 1) , 



(50) 
(51) 



where we introduce the semicolon to denote the use of the label A now. Thus, C5 for 
*(°)(E, (3; Ai) is zero and C5 for *(o)(E, (3; A2) is 



C5(S,/3;A2) = 6C5(A,/?). 



(52) 



2. The second-order energy 

To be convenient, we use the expression ^'(°)(E,/3; A3) instead of ^'(°)(E,/3, —1). For the 
state ^^'^)(E, P; Aj), the second-order energy correction can be written in the form 

C2„(S,/?;A0 

n>3 



y(^)(E,/3;A0 = -E 



(53) 
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where 



C2„(S,/3;A,)= D^^{e,L,i',L'), (54) 

e,L,£',L'>l 
£+L+i'+L'+2=2n 

D,,{e, L, f , L') = — E /a.(^, ^, L^, e, L')(G, + (2 - i)^^-)^^) . (55) 

LsLt 



and are 



{£, L, L„ Lt, e, L') = ^„(0, 0, 0, 0, £, L, L„ L„ f , L') (56) 
-^„(0, 0, -1, 1, L, L„ L„ f , L') - ^„(0, 0, 1, -1, L, L„ L^, f , L') 
-^7o(-l, 1, 0, 0, L, L„ Lj, f , L') + (7o(-l, 1,-1,1,^, ^, L., 1^^, L') 
+5„(-l, 1, 1, -1, £, L, L„ Li, L') - ^,,(1, -1, 0, 0, 1, L, L„ Lt, i\ L') 
+^„(1, -1, -1, 1, £, L, L„ L„ L') + ^„(1, -1, 1, -1, L, L„ L„ L') , 

/a, (£, L, L„ Lt, f , L') = r7o(0, 0, 0, 0, £, L, L„ L^, f , L') (57) 
+2^„(0, 0, -1, 1, £, L, L„ Li, L') + ^„(0, 0, 1, -1, L, L„ L*, L') 
+5„(-l, 1, 0, 0, £, L, L„ Lt, L') + ^ go{-l, 1, -1, 1, L, L„ L*, L') 

+ ^ go{-l, 1, 1, -1, L, L„ Lt, L') + go{l, -1, 0, 0, £, L, L„ L^, f , L') 

+^^„(1, -1, -1, 1, L, L„ Li, f , L') + ^ ^„(1, -1, 1, -1, £, L, L„ L^, L') , 



and 



L,L„Lt,e,L')= (58) 
^[5o(-l, 1, -1, 1, i, L, L„ Li, L') - ^„(-l, 1, 1, -1, i, L, L„ Li, f , L') 
-^,,(1, -1, -1, L, L„ Li, f , L') + ^„(1, -1, 1, -1, L, L„ Li, f , L')] . 



III. CALCULATIONS AND RESULTS 

In the center-of-mass frame, the Hamiltonian of a hehum atom can be written in the form 

where m„ is the nucleus mass, fx^ is the reduced mass between the electron and the nucleus, 
ri and r2 are the position vectors of the two electrons relative to the nucleus, and ri2 is 
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the distance between them. To calculate dispersion coefficients C„, we variationally evalu- 
ate energy eigenvalues and corresponding eigenfunctions of the initial states He(2 ^P) and 
He(2 ^P) with the correlated Hylleraas basis set 

{ri ri e-"^'-^^^y^Xi^,, r^)} , (60) 

where 3^/^;^(ri, is the coupled spherical harmonics for the two electrons forming a common 
eigenstate of and L^, and a and (3 are two nonlinear parameters, which are optimized by 
first calculating analytically the first-order derivatives of the variation energy with respect to 
these parameters and then using Newton's method to find their roots. The basis set includes 
all combinations of i > £i, j > £2, and A; > with i + j + k < Q, where Q is an integer 
controlling the size of the basis set. Then we generate the spectra of the intermediate states 
for the S, P, D, F, and G symmetries by diagonalizing directly the Hamiltonian of the 
helium atom in chosen basis sets. In addition, we need to transform the 2^-pole transition 
operator T^(cr), defined in Eq. (I3l), into the center-of-mass coordinates and calculate the 
reduced matrix elements of dipole, quadrupole, octupole, and hexapole transition operators 
(for details, refer to Ref. [3]). 

Table [I shows the convergence study of the nonrelativistic energy of the °°He(2 ^P) state 



with the increase of the size of basis set. Compared to the value of Drake et al. [20|], our 
result is accurate to about 15 digits. Table HTl gives the convergence pattern of C6(A,0), 
C6(n,+,0), and C6(S,0;Ai) for °°He(2 3F)-°°He(2 ^F) as the sizes of basis sets, including 
the initial state and the three intermediate states, increase progressively. 

Table Hn] presents the long-range interaction coefficients C5, Ce, Cg, and Cio for the 
He(2 iF)-He(2 system. Table CV] shows C5, Cq, Cg, and for the He(2 ip)-He(2 ^F) 
system. Table |V] lists C5, Ce, Cg, and Cio for the He(2 3F)-He(2 ^P) system. We note that, 
for the three He(n F)-He(n' F) systems, C5(n, — C^{Jl,f3; Ai), andC5(S,/?; A3) are zero, 
(75(11, +,/3) and Cg are positive, and both C^{/S.,I3) and C5(S,/3;A2) are negative. For the 
\E'^°^(S,/3; A3) states, Cg and Cio are negative in the three tables. 

Ovsyannikov obtained expressions for the Cg (dispersion) coefficients between two ex- 
cited atoms [qI and evaluated the He(2 '^P2)^He(2 '^^2) coefficients using the atomic dynamic 
polarizability obtained with a model potential for the He(2 ^P) atom. The diagonal elements 
of the long-range interaction were given in the jj representation; after transformation to the 



LS representation 



13| the results of Ref. [9| can be compared with ours. The one signifi- 
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cant discrepancy is between our results for the summed CqCE, 0; A) and that of Ref. i. In 
table |Vl] we compare the transformed coefficients of Ref. |9| with the present calculations. 

To the best of our knowledge, there are no other published results for the dispersion 
coefficients for interaction between He(2 ^P)-He(2 ^P) and He(2 ip)-He(2 ^P). 
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TABLE I: Convergence study for the nonrelativistic energy 
of He(2'^P) with the infinite nuclear mass. denotes the 
number of terms in the basis set. Units are atomic units. 



n N E{n) 

12 910 -2.133164190 779194 

13 1120 -2.133 164 190 779 246 

14 1360 -2.133 164 190 779 279 

15 1632 -2.133 164 190 779 281 2 

16 1892 -2.133 164 190 779 282 6 
Drake [2^ -2.133 164 190 779 283 202(5) 



TABLE II: Convergence study of C6(A,0), C6(n,+,0), 
and C6(S,0;Ai), in atomic units, for °°Re{2^P)-°°Ile{2^P). 
A2 3p, A35, A(pp)3p, and Aa/j denote respectively the sizes of 
bases for the initial state and the three intermediate states 
of symmetries ^5", (pp) ^P, and ^D. 

N^sp Nss N^pppp Nsp Cg{A,0) g6(n,+,0) C6(S,0;Ai) 

1360 560 1230 853 4083.026 993 80 2324.937 567 24 7296.821324 76 

1632 680 1430 1071 4083.026 997 99 2324.937 572 24 7296.821330 48 

1938 816 1650 1323 4083.026 998 66 2324.937 572 91 7296.82133149 
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TABLE III: The long-range interaction coefficients, in atomic 
units, for the He(2 ^P)-He(2 ^P) system . 





°°He(2 ip)-°°He(2 ^P) 


4He(2 ip)-4He(2 ^P) 


3He(2 ip)-3He(2 ^P) 


C5(A,0) 


-227.26160604(2) 


-227.31888382(1) 


-227.33762037(2) 


C5(n,+,o) 


909.04642409(1) 


909.27553527(2) 


909.35048146(3) 


C5(n,-,o) 











C5(S,0;Ai) 











C5(S,0;A2) 


-1363.56963613(1) 


-1363.91330289(1) 


-1364.02572224(5) 


C5(S,0;A3) 











C6(A,0) 


5578.63(2) 


5584.99(1) 


5587.07(1) 


C6(n,+,o) 


3606.63(2) 


3610.56(2) 


3611.84(1) 


C6(n,-,o) 


10831.74(2) 


10844.25(2) 


10848.35(2) 


C6(S,0;Ai) 


9682.78(1) 


9693.98(2) 


9697.63(2) 


C6(S,0;A2) 


1025.14(2) 


1032.23(2) 


1034.58(3) 


C6(S,0;A3) 


1966.062(1) 


1968.141(1) 


1968.821(1) 


C8(A,0) 


1556.9(6) 


1458.9(7) 


1426.6(6) 


C8(n,+,o) 


-170489.4(5) 


-170548.7(5) 


-170567.9(3) 


C8(n,-,o) 


1315208(3) 


1316393(2) 


1316781(2) 


C8(S,0;Ai) 


1106570(2) 


1107553(2) 


1107874(2) 


C8(S,0;A2) 


1788288(2) 


1790861(2) 


1791704(2) 


C8(S,0;A3) 


-16971.98(4) 


-17030.68(3) 


-17049.91(3) 


Cio(A,0) 


3.61343(3) X 10^ 


3.61541(2) X 10^ 


3.61607(3) X 10^ 


Cio(n,+,0) 


5.163864(4) x 10^ 


5.168052(3) X 10^ 


5.169423(3) x 10^ 


Cio(n,-,o) 


1.7334869(4) x 10^ 


1.7349273(2) x 10^ 


1.7353990(4) x 10^ 


Cio(S,0;Ai) 


2.0453997(5) x 10« 


2.0471586(3) x 10« 


2.0477345(5) x 10^ 


Cio(S,0;A2) 


7.6324701(1) X 10^ 


7.6394035(3) x 10^ 


7.6416728(2) x 10^ 


Cio(S,0;A3) 


-2.992224(2) x 10^ 


-2.996113(3) X 10^ 


-2.997383(2) x 10^ 
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TABLE IV: The long-range interaction coefficients, in atomic 
units, for the He(2 ^P)-He(2 ^P) system . 



Mass 


°°He(2 ^P)-°°He(2 ^P) ^He(2 ip)-^He(2 ^P) 




C5(A,±) 


-189.519232605(3) 


-189.513470009(3) 


-189.511579458(2) 


C5(n,+,±) 


758.07693042(3) 


758.05388004(2) 


758.04631785(3) 


C5(n,-,±) 











C5(E,±;Ai) 











C5(S,±;A2) 


-1137.11539564(3) 


-1137.08082006(2) 


-1137.06947676(2) 


C5(S,±;A3) 











C6(A,±) 


5186.93(2) 


5190.06(3) 


5191.91(2) 


C6(n,+,±) 


2990.19(1) 


2992.03(2) 


2992.640(2) 


Ce{U,-,±) 


10353.76(2) 


10359.97(5) 


10362.05(2) 


C6(S,±;Ai) 


9245.11(2) 


9250.65(3) 


9252.52(3) 


C6(S,±;A2) 


6162.84(4) 


6165.27(3) 


6166.13(4) 


C6(S,±;A3) 


1505.135(1) 


1506.081(1) 


1506.390(1) 


C8(A,±) 


-41757.6(6) 


-41832.7(4) 


-41856.3(4) 




-88728.8(2) 


-88801.9(1) 


-88825.1(2) 


C8(n,-,±) 


993884(3) 


994253(2) 


994375(2) 


C8(S,±;Ai) 


843530(3) 


843832(2) 


843931(2) 


C8(S,±;A2) 


2139396(4) 


2140392(3) 


2140712(2) 


C8(S,±;A3) 


-11571.12(3) 


-11605.17(2) 


-11616.32(3) 


Cio(A,±) 


2.32333(4) x 10^ 


2.32313(3) X 10^ 


2.32312(1) X 10^ 


Cio(n,+,±) 


5.236817(3) x 10^ 


5.238065(4) x 10^ 


5.238463(1) x 10^ 


Cio(n,-,±) 


1.224688(2) x 10^ 


1.224994(2) x 10^ 


1.225094(1) X 10^ 


Cio(S) =t; Ai) 


1.435423(2) x 10^ 


1.435806(2) X 10« 


1.435933(1) X 10« 


Cio(S, ±; A2) 


5.737754(2) x 10^ 


5.739452(2) x 10^ 


5.740007(2) X 10^ 


Cio(s, =t; A3) 


-2.018823(3) X 10^ 


-2.020165(2) X 10^ 


-2.020604(2) X 10^ 
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TABLE V: The long-range interaction coefficients, in atomic 
units, for the He(2 3p)-He(2 ^P) system . 



Mass 


°°He(2 3p)-°°He(2 ^P) 


4He(2 3p)-4He(2 3p) 


3He(2 3p)-3He(2 ^P) 


(75(A,0) 


-158.0449076079(7) 


-157.995476 2726(5) 


-157.979302 7289(2 


C5(n,+,0) 


632.179 630432(3) 


631.981905 089(1) 


631.917 210 917(2) 


C5(n,-,o) 











C5(S,0;Ai) 











C5(S,0;A2) 


-948.269445647(4) 


-947.972857636(3) 


-947.875816375(3) 


C5(S,0;A3) 











C6(A,0) 


4083.026 998 9(3) 


4082.935 140 9(1) 


4082.904 962 9(7) 


C6(n,+,o) 


2324.937 573 2(4) 


2325.200 148 8(4) 


2325.285 984 8(3) 


C6(n,-,o) 


8172.742842(3) 


8172.314300(1) 


8172.173845(3) 


C6(S,0;Ai) 


7296.821 336(5) 


7296.447426(2) 


7296.324 875(5) 


C6(S,0;A2) 


8367.590 5(5) 


8361.204 8(2) 


8359.115 7(4) 


C6(S,0;A3) 


1159.124 365 8(5) 


1159.378 065 7(5) 


1159.461 0336(4) 


C8(A,0) 


-30170.71(6) 


-30168.72(8) 


-30168.10(6) 


C8(n,+,o) 


17083.8(2) 


17044.4(2) 


17031.5(3) 


C8(n,-,o) 


726869.3(2) 


726728.9(2) 


726682.8(1) 


C8(S,0;Ai) 


618334.9(1) 


618201.6(2) 


618157.8(1) 


C8(S,0;A2) 


1835912.3(2) 


1835108.6(4) 


1834846(1) 


C8(S,0;A3) 


-8209.2(2) 


-8229.3(3) 


-8236.16(8) 


Cio(A,0) 


1537493(3) 


1536725(3) 


1536474(3) 


Cio(n,+,0) 


3.656172(2) x 10^ 


3.653559(2) x 10^ 


3.652704(2) x 10^ 


Cio(n,-,o) 


8.4139100(8) X 10^ 


8.410605(1) X 10^ 


8.409524(2) x 10^ 


Cio(S,0;Ai) 


9.846488(2) x 10^ 


9.842827(2) x 10^ 


9.841629(2) x 10^ 


Cio(S,0;A2) 


3.9464965(4) x 10^ 


3.94483100(3) x 10*^ 


3.9442865(4) x 10^ 


Cio(S,0;A3) 


-1.378107(3) X 10^ 


-1.378199(2) X 10^ 


-1.378225(5) x 10^ 
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TABLE VI: Comparison of the present results with the avail- 

n 

able results of Ref. [9.] for the Cg coefficients of He(2^P2)- 
He(2^P2) with infinite nuclear mass. The relation between 
the present results (LHS) and the results from Ref. Q] (RHS) 
are given in the first column. The symbols for the dispersion 
coefficients on the RHS represent values Ce{MA, Mb) for two 
atoms in the total angular momentum states Ja = Jb = 
and quantization axis along the inter-nuclear axis. 



Terms 


Present Ref. [9] 


C6(A,0) =C6(2,2) 


4083 4056 


i [CeiU, +, 0) + C6(n, -, 0)] = 2C6(1, 2) - 


C6(2,2) 4 666 4645 


i ELi C6(S, 0; A,) = 4C76(0, 2) + 1^(1, 1) 


- f C6(l,2) + |C6(2,2) 5 608 7026 
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